import matplotlib.pyplot as plt
import sympy as sp
import numpy as np

q = [3, 5]
t, i = sp.symbols('t,i')
eq = t ** 2 / (q[1] - (t / q[0]) ** 2)
x = np.linspace(0, 100, 10000)
y = x ** 2
plt.plot(x, y, label=r"$" + f"{sp.latex(eq)}" + r"$")
plt.legend()
plt.show()
